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Using the pruned-enriched Rosenbluth Monte Carlo algorithm, the scattering functions of semi- 
flexible macromolecules in dilute solution under good solvent conditions are estimated both in d = 2 
and d — 3 dimensions, considering also the effect of stretching forces. Using self-avoiding walks of 
up to N = 25600 steps on the square and simple cubic lattices, variable chain stiffness is modeled 
by introducing an energy penalty tb for chain bending; varying qb — exp(— eb/ksT) from qb = 1 
(completely flexible chains) to qb = 0.005, the persistence length can be varied over two orders 
of magnitude. For unstretched semiflexible chains we test the applicability of the Kratky-Porod 
worm-like chain model to describe the scattering function, and discuss methods for extracting per- 
sistence length estimates from scattering. While in d = 2 the direct crossover from rod-like chains 
to self-avoiding walks invalidates the Kratky-Porod description, it holds in d — 3 for stiff chains 
if the number of Kuhn segments tik does not exceed a limiting value n* K (which depends on the 
persistence length). For stretched chains, the Pincus blob size enters as a further characteristic 
length scale. The anisotropy of the scattering is well described by the modified Debye function, if 
the actual observed chain extension (X) (end-to-end distance in the direction of the force) as well 
as the corresponding longitudinal and transverse linear dimensions (X 2 ) — (X) 2 , (Rg x) are used. 



I. INTRODUCTION 

Small angle (neutron) scattering from polymers in di- 
lute solution is the method of choice to obtain a com- 
plete picture of the conformations of long flexible or 
semiflexible macromolecules, from the length scale of 
the monomeric units to the gyration radius of the chain 
molecules [H-[3|- Classical experiments have shown that 
the gyration radius R g of long flexible chains in dense 
melts (and also in dilute solution under Theta condi- 
tions) scales with chain length N according to the classi- 
cal random walk picture [a-Q]) Rg oc iV 1 / 2 , while in di- 
lute solution under good solvent conditions flexible poly- 
mers form swollen coils [l|-0, HI, R g oc N", with (H 
v sa 0.588 (in d = 3 dimensions) or v = 3/4 (d = 2). 
ft was also shown that on intermediate length scales 
the dependence of the scattering intensity S(q) on wave 
number q reflects the self-similar fractal structure of the 
polymer Q, S(q) cx q~ x l v (under good solvent condi- 
tions) or S(q) oc q~ 2 (Theta conditions or melts, re- 
spectively). Also the crossovers between these regimes 
when cither the temperature distance from the Theta 
point or the concentration of the solution are varied have 
been investigated and the length scales ruling these 
crossovers (i.e., diameter of "thermal blobs" or "concen- 
tration blobs" , respectively) have been identified [l|-|4( . 

The behavior gets more complicated, however, when 
chain stiffness plays a prominent role: only when chain 
stiffness is essentially due to chain thickness, i.e. when 
the effective persistence length £ p scales proportional to 
the local chain diameter D the problem can still be re- 
duced to a rescaled self-avoiding walk problem (Tol . ITl| . 
However, when £ p 3> D, one finds (in d — 3) a dou- 
ble crossover, since then short chains behave like rigid 



rods (i.e., R g oc N as long as R g < £ p ), and then 
a crossover to Gaussian random- walk like coils occurs, 
while for R g ss Rq (i?o will be discussed below) a second 
crossover to swollen coils (R g oc N u ) starts. Also this 
behavior has been established in beautiful classical ex- 
periments [HJ , as well as in recent simulations [n], [III ■ 
The standard model for semiflexible worm- like chains, 
the Kratky-Porod model [HI , can only describe the first 
crossover (from rods to Gaussian coils) but fails to ac- 
count for the second crossover to swollen coils, due to 
its complete neglect of excluded volume effects, ft turns 
out that this second crossover still is incompletely under- 
stood: while in the early experiments [12| it was sn g~ 
gested that this crossover occurs for nx = n* K with 
n* K = 50 Kuhn segments, independent of the persistence 
length, a Flory-type argument [HI, [l6| suggests that the 
crossover occurs for a polymer radius i?o oc l 2 p / D (corre- 
sponding to n* K oc (£ P /D) 2 ), while the simulations rather 
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can be described as an equivalent freely jointed chain of 
7ik Kuhn segments of length 2£ p each 
of an apparent conflict of the estimate \Y1 
with the Yamakawa-Stockmayer-Shimata theory [IB - [l9| 
for worm- like chains the result n%- = 50 was considered as 
a very fundamental problem 2^, 21 1. However, Tsuboi et 
al 2l| confirmed this estimate for another stiff polymer 
and concluded that the result n* K = 50 is compatible 
with the theory. We also note that in d — 2 there oc- 
curs a single crossover from rods to self-avoiding walks, 
any regime of Gaussian-like behavior is completely ab- 
sent |13l. [231. We emphasize however, that the results 
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of [TJliEliE O imply that a universal number 
dependent of £ p ) up to which the Kratky-Porod model 
holds in d = 3 does not exist. 
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It then is interesting to ask how chain stiffness shows 
up on intermediate length scales, that can be probed via 
the scattering function S(q). If one disregards excluded 
volume [l[ and bases the treatment on the Kratky-Porod 
model, one can show [23[ that the rod-like behavior at 
large q leads to a scattering law proportional to i.e. 
(in d — 3 dimensions) 

N£ b qS(q) = 7r+^(q£ p )- x , N -> oo , (1) 

where we have assumed that the chain has a contour 
length L = Nib, where £b is the bond length while the 
persistence length [6| is l v = (N/riK:)lb/2- While for 
Gaussian "phantom chains" (i.e., excluded volume in- 
teractions are completely neglected) the structure factor 
S(q) is readily found pHg in terms of the Debye function, 
and the only length that enters is the gyration radius 
(Rg) 1 / 2 , choosing a normalization where S(q = 0) = 1 

S(q) = 2[exp(-C) - 1 + C]/C 2 , C = q 2 (R 2 g ) = Q 2 N£ 2 b /6 , 

(2) 

for semiflcxible polymers the calculation of S(q) for 
chains with N finite is a formidable problem |24| - |39| . 
even in the absence of excluded volume effects. However, 
including excluded volume effects in the description of 
scattering of semiflexible chains is even more an unsolved 
problem: existing phenomenological approaches require 
the adjustment of many empirical parameters [37| . It 
will be one of the tasks that will be addressed in the 
present paper, to investigate S(q) for semiflexible chains 
numerically in the presence of excluded volume interac- 
tions between the effective monomers, varying l p over a 
wide range. 

In recent years also the behavior of macromolecules 
under the influence of stretching forces has found enor- 
mous interest (e.g. [iol - l55j ). in particular for the study of 
bio-macromolecules. Experimentally this can be realized 
e.g. by pulling at one end of a chain, that is anchored 
at a substrate with the other chain end, by the tip of 
an atomic force microscope [H, [H, HI, El, HI HI] but 
it is also conceivable to stretch polymers by the forces 
occurring when a polymer solution is exposed to strong 



shear flow [56l - l58j or elongational flow (59|. Of course, it 
is not obvious that it will be possible to carry out scat- 
tering experiments on such stretched chains and measure 
the structure factor (which then is anisotropic and has 
two relevant parts S\ i (qt i ), S± (q± ) since the direction of 
the scattering vector q relative to the stretch direction, 
either parallel, qu, or perpendicular, q±, matters). But 
nevertheless a theoretical investigation of S\ \(q\\), Sj_(q±) 
is worthwhile, because it gives detailed insight into the 
local structure of stretched chains, including also chains 
under cylindrical confinement |60|-|62j and this may help 
to understand problems such as transport of semiflcxi- 
ble polymers throu gh porous materials, or channels in 
nanofluidic devices [6(|, etc. Thus we shall also investi- 
gate the structure factor of stretched semiflexible chains, 
extending previous work on flexible chains [H, |44| . 

The outline of our paper is as follows: in Sec. II, we 
give a summary of the theoretical background, and in 
Sec. Ill we define our model and briefly recall the simu- 
lation methodology. In Sec. IV we present our results for 
the structure factor S(q) of semiflexible chains, for both 
d = 2 and d = 3 dimensions, in the absence of stretching 
forces. Sec. V describes the modifications of the struc- 
ture factor due to stretching, while Sec. VI summarizes 
our conclusions. The calculation of the scattering func- 
tion of random walk chains under constant pulling forces 
can be carried out analytically and is presented in an 
appendix. 



II. THEORETICAL BACKGROUND 

A. Definitions 

We consider here the scattering from a single polymer 
chain, assuming that the chain can be described by a 
sequence of N + 1 (effective) monomers at positions rj , 
j = 1,2,..., N + 1, so that we can define N bond vectors 
Sj = fj + i — rj. In the absence of stretching forces, the 
structure factor S{q) does not depend on the direction of 
the scattering vector q, and can be defined as 
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Note that we have chosen here a normalization for which one chain end in the +x-direction, the structure factor 
S(q -+ 0) = 1. When a stretching force is applied to becomes anisotropic. In d = 3 dimensions, the confer- 
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mations of chains still have axis-symmetric geometries, where q is oriented perpendicular to it. So we define 
and we must distinguish between 5*11(911), where q is ori- fj — (xj , yj ,Zj) — (xj , pj ) to obtain 
ented in the x-direction parallel to the force, and S±(q±), 
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Note that in d = 2 dimensions we have fj = (xj , yj ) and 
then q± ■ pj in Eq. (O needs to be replaced simply by 
lA-Uji °f course. 

We also stress that in this paper we are not at all con- 
cerned with effects due to the local structure of (effec- 
tive) monomers, such as e.g. chemical side groups, etc.; 
such effects show up at large q when one considers the 
scattering from real chains [29(. We next define our no- 
tation for characteristic lengths of the chain. Assum- 
ing a rigidly fixed bond length £{, between neighboring 
monomers along the chain, the contour length L is 



N£ h 



(6) 



The mean square end-to-end distance (in the absence of 
stretching forces) simply is 



(R 2 




(7) 



with dj = fj + i — fj being the j-th bond vector, and the 
mean square gyration radius is given by 
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where fcM = X) l{N + 1) is the center of mass posi- 
3=1 

tion of the polymer. 

In the presence of stretching forces, the chain takes 
a mean extension (X) and mean square extensions also 
become anisotropic, 
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Eq. (fTU)) refers to the three-dimensional case, for d — 
2 the second term in the right hand side needs to be 
omitted. A related anisotropy can then be stated for the 
gyration radius square as well, namely 
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for d = 3, again the term (zj — z^) 2 simply is omitted in 
the case d = 2. 

We also recall that the mean square gyration radii de- 
scribe the scattering functions at small q. So in the ab- 
sence of stretching forces we have 

S(q) = l-q 2 (R 2 g )/d. / q^O, (13) 

while if stretching forces are present, one finds instead 

S||te||)=l-gfl<i£n>, q^Q, (14) 



S x ( qi _) = l-ql(Rl x )/(d-l) 



0. 



(15) 



In addition to the limit q — > 0, also the limiting behavior 
of S(q — > 00) is trivially known: then all interference 
terms in Eq. ([3]) average to zero, and only the terms 
j = k contribute to the double sum. Hence we obtain 



S(q — > 00) = 



1 



N + l' 



S\\(q\\ -+ 00) = Sj_(qx^>- 00) 
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(16) 
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irrespective of the value of the persistence length £ p , the 
value of an applied force /, etc. 

At this point we emphasize, however, that it is not 
possible to write down a general definition for the per- 
sistence length £ p that would be both universally valid 
and practically useful Hoi [ill [l3| . The definition often 
found in textbooks [H, l63j in terms of the asymptotic de- 
cay of the bond vector orientational function (of a very 
long chain, N — > oo) 

(cos 0(s)) = (dj ■ dj +s )/(3j ■ dj) oc exp(— s£b/£ p ) , s — >• oo 

(17) 

makes sense only for Gaussian PHANTOM chains, and 
is not applicable to real polymers under ANY CIRCUM- 
STANCES, since the asymptotic decay of (cos0(s)) with 
the "chemical distance" s£b along the chain always is 
a power-law decay. In fact, for £ p <C s£b <C L one 
has [T3, [ill [TE E3467I 



(cos 0(s)) oc s 13 



(18) 



where /? = 3/2 both in melts [64|, [65| and for chains in 
dilute solution under Theta conditions [l(| [HI . For the 
case of good solvent conditions, which is the problem of 
interest for the present paper, one rather finds the scaling 
law [H M 



= 2(1-1/) 



(19) 



which yields (3 = 0.825 (d = 3) and = 1/2 (d = 2), 
respectively. In simple cases, such as the semifiexible 
extension of the self- avoiding walk model studied in 13 1 
and further investigated in the present work, one can 
rather use an analog of Eq. (fT?) but for short chemical 
distances, 



in d = 3 dimensions [13j. Unfortunately, Eq. ([2"Tj) is not 
straightforwardly applicable for chemically realistic mod- 
els (such as alkane chains when £b means a bond between 
two successive carbon atoms, but the all-trans state cor- 
responds to a zig-zag configuration with a nonzero bond 
angle 0(1)). It also is not useful for coarse-grained models 
of polymers with c omp lex architecture, such as bottle- 
brush polymers [ToL lll| . Thus we emphasize that for our 
model Eq. (f2~Tj) is a practically useful definition, while 
for real polymers studied experimentally the estimation 
of £ p is a delicate problem. The same caveat applies 
for the Kuhn length £r, which is £k = 2£ P for worm- 
like chains, but the latter does not apply in solutions, 
as stated above. In dense melts, £x/£b — 6(R 2 )/(£ 2 N) 
is supposed to hold, but due to local interactions with 
neighboring monomers in a dense environment it is not 
obvious that for a melt is a relevant parameter for a 
chain under good solvent conditions. 

Thus, it is a clear advantage of our model calculations 
that via Eqs. (|2U|) . (j2"TT) accurate direct estimates of £ p 
are possible, unlike in experiment. These estimates for 
£ v are tabulated in Ref. fl3i. 



B. Theoretical Predictions for the Scattering 
Function of Single Polymers in Good Solvents in the 
Absence of Stretching Forces 

The classical result for the scattering from Gaussian 
chains is the well-known Dcbyc function [lHj] 



SbcbycC?) 



l - 



l 



q 2 (R% 



■[l-cxp(- g 2 (i? 2 g ))] 



Note that for large q this reduces to Sbcbye^) 



(22) 



(cos 6(s)) = exp(-s£ b /l p ) , < s < £ p /£ b (20) 2/(q 2 (B*}), reflecting the random-walk like fractal struc- 



Eq. ([2TJ]) is useful for the simple model that will be stud- 
ied in the present paper, namely the self-avoiding walk 
(SAW) on square and simple cubic lattices with an energy 
penalty e& for "bond bending" (i.e., kinks of the SAW by 
90 degrees) , which then serves as a convenient parameter 
to control the persistence length. Since for eb/ksT < 2 
the power law, Eq. flTBl . already starts to set in even 
for small s of order unity already, we use in practice an 
alternative definition, 



ture of a Gaussian coil, SD e b ye (q) oc q~ x l v MF with Umf 
1/2. Eq. (|22p does not tell how large q can be in order 
for this power law to be still observable. For semifiex- 
ible Gaussian chains, the contour length L = N£b can 
be written as L = n p £ p and the mean square end-to-end 
distance and gyration radius are 



(R 2 ) 



-[1 - exp(-n p )] , 



(23) 



V4 = -1/M<cos0(l))), 



(21) 



which is equivalent to Eq. (|2U|) . if Eq. (|2"0"|) holds over a 
more extended range of s. An alternative method con- 
siders the distribution function P(n str ) of n s t r succes- 
sive bond vectors a$ having the same orientation without 
any kink, which is found to behave as [13( P{n s tr) — 



a p exp(~n str /n p ), with 



h p , I V p 



being constants. Both 



n p can be taken as alternative estimates of a 



{n str ) and 

persistence length, and for large n p these estimates agree 
with the result from Eq. (|2~T1) to within a relative accu- 
racy of a few percent (or better), both in d = 2 and 



6(R 



i 
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Tip Tip Tip 



[1 - exp(-n p )] (24) 



From Eqs. (|2"5|) . one can clearly recognize that Gaus- 
sian behavior of the radii is only seen if the number n p 
of persistence lengths that fit to a given contour length 
of the chain is large, n p 3> 1 (for n p of order unity, a 
crossover to rod- like behavior occurs). Since g _1 also is 
a length scale, one concludes that the Gaussian coil be- 
havior reflected in Eq. (|22|) also implies that a scale g _1 
requires that a subchain with this gyration radius con- 
tains many persistence lengths as well, i.e. Eq. (|2"2"|) can 
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only hold for 



which for large q varies like 



In the regime 



q£ p -C 1 . 



t- 1 « q « t; 1 



(25) 



(26) 



we expect that the scattering function will resemble the 
scattering function of a rigid rod of length L vo d [69| 



S m d(q) 



qL IO d 



qL Ta d 



dx 



sin a; 1 — cos(qL IQ d 



qL m d 



(27) 



SVod(<? -> OO) = 7r/ (gLrod) • 



(28) 



Note that Eqs. (|27). (|28|) refer to a rigid rod on which 
the scattering centers are uniformly and continuously dis- 
tributed. In the lattice model that is studied here, the 
scattering centers are just the subsequent lattice sites 
along the rod, and since for a rod of length L TO d there 
are then L ro d + 1 such scattering centers, one has [70l ] 



irod 



S m d(q) = 



L ro d 



-1 



k=0 



rod 



l - k) 



sin qk 
qk 



q < 7r . 



(29) 



Both Eqs. (|27l) . (|2^|) have a simple smooth crossover from 
Srod(g) - 1 - <7 2 (i?2) rod /3 with (i?2) rod = L 2 m J\2 to the 
l/<7 power law (Eq. P5)l ). Of course, on the lattice con- 
sideration of g > 7r does not make sense, since distances 
of the order of a lattice spacing and less are not mean- 
ingful. 

While for a Gaussian coil no direction of q is singled 
out, for a rod it makes sense to consider also the special 
case where the wave vector q is oriented along the rod; 
rather than considering the case when all orientations of 
q are averaged over, as done in Eqs. (l27l) - (j29|) . Then 
one rather obtains, <?|| being the component of q parallel 
to the axis of the rod [56| 



can be cast in the form 



(<7||£rod) : 



r[l-cos(<Z||L„d)]. (30) 



([50)1 leads to an oscillatory decay 
= 1 for q^> = fc(27r/£ rod ), k 
ami S I0 d(q\\) hence has zeros for all 

L T 



Note that Eq. 
since cos(g||L ro d] 
0,1,2, L ro d 

(k) 

9|| , fc = 0, 1, 2, . . . , iv ro d 

When we consider the scattering from semiflexible 
Gaussian chains, we now expect a smooth crossover be- 
tween the Debye function, Sbebye (<?) and the rod scatter- 
ing, Eq. (|27|). similar to the smooth crossovers from rods 
to Gaussian coils, as described for the radii by Eqs. (|2"5)l 
and (|24|) . It turns out that this is a formidable prob- 
lem, and no simple explicit formula exists, despite the 
fact that excluded volume effects still are neglected fl7i — 
IHLIIS-Ial, Kholodenko [13 derived an inter- 

polation formula which describes the two limiting cases 
of Gaussian coils and rigid rods exactly, and which is 
expected to show only small deviations from the exact 
result in the intermediate crossover regime. His result 



S(q) = -[h(x) - -I 2 (X)], 



3L/2£ P 



where 



I n (x) = J dzz n - L f(z), 



and the function f(z) is given by 



(31) 



(32) 



/(*) = 



1 sin\\.{Ez) 
E sinh z 
1 sin(Ez) 
E sinh z 



q< 
q > 



3 
2f,, 

3 
2f„ 



(33) 



with 



E = 



1 



2qi p 



1 1/2 



^ 2q£ p 



- 1 



1/2 



(34) 

In addition, Stepanow [35| has developed a system- 
atic expansion of the scattering function in terms of the 
solution for the quantum rigid rotator problem, which 
converges fast if L/£ p is not too large. Note that the op- 
posite limit, Lltn — » oo, has already been considered by 
des Cloizeaux [23[, and his result has been quoted in the 
introduction (Eq. (H}). This expansion (as well as equiv- 
alent representations written as continued fractions [33j) 
can only be evaluated numerically [7l| . 

However, a few qualitative statements can be made 
on the structure factor S(q) in the representation of a 
Kratky plot, qLS(q) plotted vs. Lq. While Eqs ([27])- 
(1291) imply a monotonous increase from the straight line 
qLS(q ~ 0) = qL towards the plateau value qLS{q ^> 
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2tt / L) = 7T for simple rigid rods, for semiflexible poly- 
mers this Kratky plot exhibits a maximum, since in the 
regime of interest we may crudely approximate S(q) by 
its leading terms 



qLS(q)nqL-hqL) 3{U; > 



V- 



(35) 



q max = l/^{R 2 ),q max LS(q max ) = ^L/^{Rj) (36) 
i.e. From Eqs. (|3"5|) . (|2"H) we immediately find, for n p 3> 1, 



<?ma X = y3/(£ p L) , q ma , x LS{q ma ^) = — y L/£ p . (37) 

Using the full Debye function one finds a different prefac- 
tor, g max « -\/QA/(lpL), but the general scaling behavior 
is the same as given by Eq. Q37p. Thus, when L is known 
(as is the case in simulations) for semiflexible Gaussian 
chains estimation of the coordinates of the maximum in 
the Kratky plot allows a straightforward estimation of 
the persistence length £ p . 

Considering now the effects of excluded volume, we em- 
phasize that Eq. (|3"6"f still is supposed to be valid, while 
Eq. (l37l) no longer holds. In particular, in Ref. Il3| it 
was shown that in d = 2 the Kratky-Porod model [lj] of 
semiflexible chains, on which Eqs. (|2"3"|) . ([Ml) , and (|3"Tj) - 
(|34|) are based, has no validity whatsoever: rather it was 



shown that around n p — 1 a smooth crossover from the 
rigid rod behavior to the behavior of two-dimensional 
self-avoiding walks occurs. Thus, we expect similarly in- 
stead of Eq. (p4|) that (recall n p = L/£ p ) 



= /(«p)> witri f( n p <l) = n 



p i 



(38) 



C[/Cf and C|/C| are expected to be universal, how- 
ever (for Gaussian chains (R 2 ) / (R 2 ) = 6), for both flexi- 
ble and semiflexible chains. 

In d = 3 dimensions, however, the situation in the 
presence of excluded volume is considerably more in- 
volved [la ]. For the end-to-end distance (R 2 ) of semi- 
flexible chains, two successive crossovers were found: for 
n p sa 1 a crossover from rods to Gaussian coils occur, 
while excluded volume effects become prominent for 



n„ > n r 



n* p ex (£ P /D)i 



(43) 



where we have introduced the chain diameter D as an- 
other characteristic length that may be needed in general 
(while in our model D — £b, however), and £ is an expo- 
nent that is not yet known precisely. Arguments based 
on Flory theory yield [13j, ll5L [l6| ( = 2 while Monte Carlo 
results rather suggested [13| £ w 1.5. We recall that in 
d = 3 Flory arguments are not Q exact, implying 0, Q 
v = 3/5 instead of v = 0.588. A similar double crossover 
from rods to first Gaussian coils and then to d = 3 self- 
avoiding walks is expected to be visible in (Rg), too. If 
we could rely on Flory theory, we would predict from 
these considerations that 



(R 2 g )/(£ P L) - C* g {n v lnlf v - X , n p > n* p (44) 
and hence (using the Flory value v = 3/5) 



(R 2 .) = C s q L 6 ' 5 {i p Df' 5 



(45) 



From q n 



= (y^I/V^pD) 1 / 5 )" 1 the persistence 
length £ p can be inferred, provided Cg has been deter- 
mined. However, if £ p 3> £(, and D — £{,, one can study 
the regime 1 < n p < n*, where Gaussian statistics for 
the gyration radius is still applicable, and hence Eq. (|37|) 
applies. 



as in Eq. (J24J) , but 



f(n p » 1) = C, 



, ^2 = 3/4, 



(39) 



Cg being (for £ p /£b ~> 1) a non-universal constant. 
Eqs. (13"B"1) . (HHJ) hence imply in d = 2 dimensions 

= C s g ll /2 L 3/2 , L^oo, (40) 
in full analogy to the result for fully flexible chains, 

{R 2 a )=Cf£l /2 L^ 2 , L^oo, (41) 

where is another non- universal constant. Similar re- 
lations were found [13( for the end-to-end distance 

(R 2 ) = Cl£l /2 L 3 / 2 (flexible) , 

(R 2 ) = C s Jl /2 L 3 / 2 (semiflexible) (42) 

with C[ and C| other (non-universal) constants for flex- 
ible and semiflexible chains, respectively. The ratios 



C. The structure factor in the presence of 
stretching forces 

For Gaussian chains under stretch, where a force / 
is applied at a chain at one end in the +x-direction, 
the other end being fixed at the coordinate origin, the 
structure factor S(q) has been derived by Benoit et 
al. [40| as follows {q = (q x ,qy,qz) with <jfii = q x and 



i 2 A 



S l{ ( qil ) = 2Re 



S±(q±) = 2 



exp(-Xi) - 1 + X 



II 



exp(-Xx) -l + X ± 



(46) 



(47) 



where X x = q]_{{R 2 ) / QX^) , with (R 2 } the mean- 
square end-to-end distance of the chain in the absence 
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of any force (/ = 0), and Aj_ describes the modification 
of the Gaussian distribution in the transverse directions 
(y and ^-direction, for d = 3). The quantity Xu is com- 
plex (therefore the real part of Eq. PS|) is taken) and is 
given by 



sian chains) that A^ = 

(i«,»/(«>o + (R* 



(3/2)(iJ» )/(i? 2 } 



expected. 



m v ) + 

Similarly, 



Eq. (gBJ) yields for q^ ->• 0, using also Eq. (fT4"| 



ATii 



<Z 2 |(i? 2 )o 
6A" 2 



(48) 



where An describes the modification of the Gaussian dis- 
tribution in the x-direction (parallel to the force). Benoit 
et al. [loj explicitly state that their result is restricted 
to deformations of small amplitudes, and do not specify 
how A_l 5 Ai are related to the applied force. However, 
considering the small q expansion of Eqs. (l46l) . (|47)) one 
can relate these parameters to the mean square gyration 
radius components of the chain, since for X± -c 1 



S(q±) = 1-«1 



2 (R 2 ) 



18A2 2 



= i-q 2 ±(R 2 g i>/2, (49) 



where in the last step Eq. (fT5")) was used. Hence 
we conclude (note that (Rg)o — (R 2 )o/6 for Gaus- 



S(q\\) = 1 



(R 2 



i8Ar 



12 



9,\U 



(50) 



and hence we see that An in X\\ can be expressed in 
terms of the gyration radius component (R 2 ^ of the 
stretched chain and the extension (X). However, since 
from the work of Benoit et al. [40] it is not clear that 
Eqs. (146[) . (|4"T|) are applicable for conditions where (X)/L 
is not very small, we hence rederived Eqs. (|46|) . (l47|) by 
an independent method, which is more transparent with 
respect to the basic assumptions that are made. This 
derivation is presented in an Appendix, and it shows that 
S(q\\) can be cast into the form 



%||) 



(N + iy 



exp 



I- J 

N 



cos 



\i-3\{X) 
N 



(51) 



r 



which is equivalent to Eq. (|46[) but with a somewhat dif- 
ferent expression for X\ i , namely 



X 



1 



|| = -^((x 2 )-(X) 2 )+ i(? ||(X>. 



(52) 



It is interesting to note that Eqs. (|46|) and ([STp can be 
given a very simple physical interpretation: with respect 
to the correlation in stretching direction, the stretched 
polymers is equivalent to a harmonic one-dimensional 
"crystal" (which at nonzero temperature lacks long range 
order, of course) of length Na = (X), a being the "lattice 
spacing" of the crystal. 

Writ ing the Hamiltonian of the one dimensional chain 



as 



*4£ 



El 

m 



{x e+ i - xi - a) 2 



(53) 



where point particles of mass m have positions X£ and 
conjugate momenta 71^, and the spring potential coupling 
neighboring particles is written in terms of the sound ve- 
locity c. At T — 0, particles would be localized at po- 



sitions Xn^ 



,<°) 



+ na, n 



0,1, 



N. So it makes 



sense to consider displacements relative to the ground 



state, u rt 



(0) 

putting x a 



at the 



origin. Due to the harmonic character of this "crys- 
tal" , one can calculate the mean square displacements 
easily to find (assuming periodic boundary conditions) 
that ({u n — wo) 2 ) = na 2 ksT/ (mc 2 ) = nS 2 , where 5 char- 
acterizes the local displacement for two neighboring par- 
ticles. Applying the formula also for the end-to-end dis- 
tance of a chain without periodic boundary conditions, 
((un — uo) 2 ) = NS 2 , one immediately finds that S(q\\) 
for the harmonic chain yields the above expressions of 
S(q\\), since X = xn — xq — Na + ujy — Uo, (X) = Na, 
and (X 2 ) = N 2 a 2 + NS 2 = (X) 2 + N5 2 . This considera- 
tion also emphasizes that a condition (X) <C L = Nif, in 
fact is not required for the validity of Eqs. (141)1) - (| 5 2 1) . 

For the unstretched case {(X) = 0) Eqs. (pg ]l -([5"2" ]l re- 
duce to Eq. (|2"2"|) . as it should be. We recall that accord- 
ing to the Kratky-Porod model simple approximations 
for the extension (X) of a chain as a function of the force 
can be derived (see [1 3j j for a review), namely 



ftp 


3 (X) 


1 


k B T 


4 L 


8(l-(X)/L)2 


and 






ftp 


3 (X) 


1 



1 



k R T 4 L 



4(1 - (X)/Lf 



1 

4 



d=2, (54) 



d = 3. (55) 
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TABLE I: Values of persistence lengths £ p /£b for semiflexible chains in d = 2 and d — 3, estimated by Eq. (|21[) . and the crossover 
length N* between the intermediate Gaussian regime and the SAW regime in d = 3, estimated empirically from Fig. 7b of 
Ref. (N* = 7V rod = 2£ p /£ b in d = 2) for various values of q b . 





0.005 0.01 0.02 0.03 0.05 0.10 0.20 0.40 1.0 


e p /h(d = 2) 
e p /e b (d = s) 

N*(d = 3) 


118.22 59.44 30.02 20.21 12.35 6.46 3.50 2.00 1.06 
51.52 26.08 13.35 9.10 5.70 3.12 1.18 1.12 0.67 
36000 9000 1850 700 180 41 11 - 



Eqs. (jM]) . (I5"5l) imply in the linear response regime, where 
(X) cx /, that 



ft P d (X) 



k R T 2 L 



(56) 



However, from linear response one can show generally 
that 



(X) = f(X 2 ) /(k B T) = f(R 2 ) /(dk B T) . 



(57) 



where (R 2 }o is the mean square end-to-end distance in 
the absence of forces. Eqs. ([56")) , ([57)1 are compatible with 
each other for Gaussian semiflexible chains, for which 
(R 2 )a = 2£ p L {Eq. (|23|)}. but are incompatible in the 
presence of excluded volume forces. In this case, one 
observes a crossover from the linear response regime, as 
described by Eq. ([57]) together with Eq. (g5J| for d = 2 
and a result analogous to Eq. (J55j, namely 

(R 2 ) = C s e L^ 5 (£ p D) 2 ^ 5 , d = 3 (using v » 3/5) , (58) 

to the so-called "Pincus blob" [4l| regime, described by 
a power law for the extension versus force relation 



(X)/L^x(f£ p /k B T) 1 ^- 



(59) 



While in d = 2 Eq. (HJj holds up to (A") /L of order 
unity, where then saturation effects ((X) / L — > 1 for 
large enough /) set in, in d = 3 the regime of validity 
of Eq. (1571) is much more restricted, namely we have to 
require [13j 



£ p = k B T/f > R cx £ 2 /D 



(60) 



For stronger forces (corresponding to £ p < i?o) the 
Kratky-Porod results, Eqs. ([54)) , ([55]) , are expected to 
become valid. In the Pincus blob regime, also nontriv- 
ial power laws for the fluctuations (X 2 ) — (X) 2 and the 
transverse linear dimensions are predicted |54j 



(X 2 ) - (X) 2 cx (Rl) cx {]l P lk B Tfl v - 2 . 



(61) 



Since we are not aware of any treatment of the structure 
factor of the Kratky-Porod model under stretch, we shall 
use Eqs. ([55]). fljTJ also for semiflexible chains (but using 
the numerical results for (X) and (A" 2 ) — (AT) 2 , (R 2 ± ), 
rather than theoretical predictions). 



III. MODEL AND SIMULATION TECHNIQUE 

Our model is the standard self-avoiding walk (SAW) on 
the square and simple cubic lattices, effective monomers 
being described by occupied lattice sites, connected by 
bonds. Each site can be taken only once, and thus we 
realize the excluded volume interaction. The lattice spac- 
ing henceforth is our unit of length, £{, = !. We intro- 
duce an energy e b for any kink the walk takes (by an 
angle of ±90°). Any such kink introduces hence a factor 
qb = exp(—€b/k B T) to the statistical weight of the walk. 

In the presence of a force / coupling to the extension 
X of the chain in x-direction, the statistical weight gets 
another factor b x , with b — exp(/ /k B T). Then the parti- 
tion function of a SAW with N bonds (i.e., N + l effective 
monomers) and A^ond local kinks becomes 



(q b ,b) = C(N,N benA ,X)c 



N bmdh X 



(62) 



config. 



By the pruned-enriched Rosenbluth method (PERM) it is 
possible to obtain estimates of the partition function and 
quantities derived from it (e.g. (X), (A 2 )) and additional 
averages such as S(q), using chain lengths up to A = 
25600. Both the chain stiffness and the force / have 
been varied over a wide range; for q b = 1 one has fully 
flexible self-avoiding random walks, while for qi, = 0.005 
the persistence length (computed from Eq. ([2T]) ) is of 
the order of 120 in d = 2 and 52 in d = 3 (Table I 
lists our corresponding estimates). For technical details 
on the im plem entation of the algorithm, we refer to the 
literature [13j. 



IV. RESULTS FOR THE SCATTERING 
FUNCTION OF UNSTRETCHED CHAINS 



We start with our data for the mean square gyration 
radius (R 2 ) , normalized by the square of the Kuhn length 
Ik = 2£ p , plotted vs. the number of Kuhn segments tik = 
L/£k = N£b/(2£ p ), Fig. la since this was the represen- 
tation chosen for the experimental data of Norisuye and 
Fujita 12] , which we reproduce in Fig. lb. Both diagrams 
show the same range of abscissa (30 < tik < 3000) and 
ordinate (5 < (R 2 g )/{2£ p ) 2 < 1000). The qualitative sim- 
ilarity between both simulation and experiment is strik- 
ing. Since only the regime of rather large tik is shown, 
the crossover from rods to Gaussian chains is not included 
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FIG. 1: (a) Log-log plot of normalized gyration radius square (R 2 ,) / (2£ p ) 2 versus the number of Kuhn segments, for the range 
30 < riK < 3000, for chains of widely varying stiffness, and comparison to corresponding experimental data (b), taken from 
Norisuye and Fujita 12]. The full straight line is the Kratky-Porod model, Eq. (I24|l . 




FIG. 2: (a) Log-log plot of (R 2 g )/{2£ P L) versus uk for d = 3, and including data for widely varying persistence lengths. Note 
that the worm-like chain result, WLC (Eq. (|24[) ). describes correctly the crossover from rods to coils, but not the onset of 
excluded volume effects, (b) Same as (a), but choosing N/N*(qb) as an abscissa variable (N*(qi,) was already estimated for a 
similar scaling plot for the mean square end-to-end distance, cf. Fig. 7b of Ref. Here L = Nti,. 



(the full straight line represents the Gaussian chain be- 
havior, as described by Eq. for tik = (l/2)n p » 1). 
One can see that Eq. (|24[) works for very stiff chains and 
not too large tik, while for large nx systematic devia- 
tions occur, which can be attributed to excluded volume 
effects. Both the simulation and the experiment include 
data for widely varying persistence lengths (in the exper- 
iment, this could only be achieved by combining data for 
chemically different polymers in this plot). From their 
results (see Fig. lb) the experimentalists concluded that 
the excluded volume effects set in for tik > 50, irrespec- 
tive of the precise value of the persistence length. 

However, this latter conclusion needs to be questioned: 
in fact, for large tik the data do not superimpose in 



this representation for different choices of £ p , indicating 
that the behavior is more complicated. To elucidate this, 
we take out the leading power law in the Gaussian coil 
regime, plotting (Rg)/(l p L) versus njc over the full range 
(Fig. 2a). One sees that nice scaling behavior occurs with 
respect to the crossover from rigid rods to Gaussian coils; 
in this regime, Eq. (|24p works in d — 3. However, now 
one can see rather clearly that the crossover from Gaus- 
sian coils to SAWs does not scale in this representation: 
rather for large tik the curves "splay out", the larger 
£ p the longer the data follow Eq. (f2~4")) . before an onset 
of excluded volume effects can be seen. This behavior 
has already been studied in Ref. [l3[ with respect to the 
end-to-end distance. Empirically, it was found that seal- 
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FIG. 3: Log-log plot of (R 2 g )/(2£ P L) versus n K for d = 2, 
including data for widely varying persistence lengths. Note 
that in d — 2 there occurs a direct crossover from rods to 
SAWs, an intermediate Gaussian regime in absent. 



ing N with N* , where N* oc £ 2 p b rather than N* = l p . 
Fig. 2b shows that a master curve results as an enve- 
lope of the curves for individual £ p . We also recall, that 
Flory arguments predict N* oc £ p , cf. Eq. (|43t and the 
subsequent discussion. The data in Fig. 2 are fully anal- 
ogous to our data on the end-to-end distance that were 
discussed recently elsewhere [l3[. In d = 2, however, the 
behavior is clearly simpler (Fig. 3): there occurs a sin- 
gle crossover from rods to d = 2 SAWs, and a regime 
where the Kratky-Porod worm-like chain model presents 
a faithful description of the data is absent. These results 
for (R 2 ) confirm our earlier analogous findings [l3[ for 

{RD- 

Fig. 4 shows some of our raw data for the structure fac- 
tor S{q). For small q, one recognizes the Guinier regime, 
S(q) « exp(-<7 2 (i?2)/3) 1 - q 2 (R 2 g )/3, and then a 
crossover occurs to the power law of SAWs or of Gaussian 
chains (the latter is seen clearly only for d — 3 and very 
stiff chains). At large q and stiff chains, the expected q~ l 
behavior is in fact compatible with the data. 

It turns out that an analysis of S{q) in the form of 
Kratky plots (Eq. ([35])) is more illuminating, cf. Fig. 5: 
The location of the maximum in the Kratky plot, as dis- 
cussed in Eqs. (|3"tJ|) . (|3"T1) . is easily identified, and it shows 
the expected scaling with L/£ p both in d — 3 and in d = 2 
(Fig. 6). In d = 3, one notes that with increasing L/£ p a 
crossover from Gaussian behavior to SAW behavior oc- 
curs. Our data also confirm that for rather stiff chains in 
d = 3 both the Kholodenko and the Stepanow theories 
describe S(q) very accurately. Note that Fig. [5p refers 
to rather short chains, for which strong effects due to 
excluded volume interactions are not yet expected, and 
hence the good agreement with the theories is not sur- 
prising. 

Another issue of interest is the behavior qLS(q) in the 
rigid rod limit, where one clearly notes the approach to 
7r (Fig. 5). Can we then use S(q) in this region, applying 
the des Cloizeaux formula {Eq. ([!])} to extract quanti- 
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FIG. 4: Structure factor S(q) plotted q, on log-log- scales, for 
d = 2 (a) and d = 3 (b); only data for L — 25600 are shown, 
but many different choices of the stiffness parameter ^i, are 
included. The straight lines indicate the rod-like behavior at 
large q (slope = — 1) and the SAW behavior for flexible chains 
(slope = —1/u). Also the slope expected in the Gaussian 
regime is included (slope = —2). Only data up to q = it have 
been included (larger q cannot be studied due to the lattice 
character of our model). 



tatively reliable estimates for the persistence length £ p 
from a plot of qLS(q) versus 1/q ? Fig. 7 suggests that 
although a regime occurs where the variation is linear in 
{q£ p )~ l , the coefficient of this linear variation is incon- 
sistent with the des Cloizeaux result. We have no final 
answer to offer to explain this discrepancy; we suspect 
that in the regime where q^ 1 and £ p are of the same or- 
der, the discreteness of our lattice model (opposed to the 
Kratky-Porod continuum model) might matter. 



V. RESULTS FOR THE SCATTERING 
FUNCTION OF CHAINS UNDER STRETCH 

While for unstretched chains it is only (R 2 g ) as a mea- 
sure of the linear dimension of the whole chain which 
is relevant (Figs. 1-3), for chains under the action of 
stretching forces anisotropy of the conformation of the 
chain comes into play. However, from the small q ex- 
pansion of Eqs. ((1]), §5§ one can show straightforwardly 
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FIG. 5: (a) Rescaled structure factor qLS(q) plotted against Lq for d = 3 and L = 25600, including 5 choices of the stiffness. 
The result for Gaussian chains (Debye function) and for continuous rigid rods (for which qLS(q) — > tv for large q. cf. Eq. (1281) ) 
are included for comparison. Also predictions obtained from the formulas proposed by Kholodenko {Eqs. (|31[) - (|34|l } are shown, 
(b) Rescaled structure factor qLS(q) plotted against q£ p for d = 3 and for L — 400, including both the predictions due to 
Stepanow [3J| and Kholodenko {Eq. (|31[) - (|34[l }. which are essentially indistinguishable on the scale of the figure, (c) Same as 
(b), but for 3 different choices of L for = 0.005. 



that S'||(g'||) yields information on (R 2 ^) and S±(q±) on 
(Rg j_) cf. Eqs. (fl4|) and (|15[) . Since for larger g|| where 
also the extension (X) of the chain along the direction 
of the force enters the description of the scattering, we 
begin with a discussion of these linear dimensions and 
describe their variation as a function of the force / (for a 
more detailed discussion and related results, we refer the 
reader to Ref. 13]). 

Fig. 8 shows typical data of (Rg ±) versus fip/ksT, 
both for rather flexible chains (qb = 0.4) and for rather 
stiff chains (qb = 0.01 in d = 2, <ft = 0.05 in d = 3, re- 
spectively). We recognize three regimes: for very small 
forces (R 2 ± ) « (R 2 ± )o, the unperturbed value in the 
absence of forces. In this linear response regime, the 
force orients the coil without deforming it. Then we rec- 
ognize a regime where (R 2 j_) decreases according to a 
power law, namely Eq. (1611) . This power law holds in the 
regime where the radius of the Pincus blob, £ p = ksT/f, 
is smaller than the unperturbed radius, but much larger 
than the persistence length l v itself. Thus this is the ana- 
log of the "Pincus blob" law that yields another power 



law for the extension vs. force curve, Eq. (1591) . The phys- 
ical picture invoked here for the chain is an elastic string 
of Pincus blobs, (R 2 ± ) describing the mean square dis- 
placement of this string in transverse directions. As one 
can see from Fig. [FJ the data indeed are compatible with 
the predicted power law in d = 2, and in d — 3 at least 
for the flexible chains. For stiff chains in d = 3, the 
regime where Eq. (f6Tj) holds is more restricted, since the 
Kratky-Porod regime has a more extended regime of va- 
lidity, effects due to Pincus blobs can only be detected in 

a regime (R 2 g ^)l /2 > £ P > R* <x £ 2 p /D, cf. Eq. (60) [ljj. 
Therefore we have not included very stiff chains in Fig.|5jD 
(for qb — 0.005, leading to i v ps 52£b, excluded volume 
effects, which also are responsible for the existence of 
Pincus blobs, could even for chains as long as N = 25600 
hardly be detected in the chain linear dimensions in the 
absence of a force, cf. Figs. 1-3). So this failure to detect 
Pincus blobs for very stiff long chains in d = 3 dimen- 
sions is hardly surprising (although Eq. (|61l) ultimatively 
will become valid as N — ¥ oo, irrespective how large £ p 
is). 
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FIG. 6: Log-log plot of the inverse position of the maximum in the Kratky plot, (^ p q ma x) 1 , versus the rescaled contour length 
L/£ p , for d — 2 (a) and d = 3 (b). The full straight line included in (a) is described by the equation 0.19(L/l p ) 3 ^ 4 '. In case 
(b), two straight lines with slope v = 1/2 and v — 0.588 are included, to illustrate the crossover. In both cases the persistence 
length tp is varied over a wide range, and the contour lengths chosen are L = 400, 1600, 6400 and 25600, as indicated. 
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FIG. 7: Plot of qLS(q) vs. (qLp)" 1 , for d = 2 (a) and d — 3 (b). All data are for L = 25600 only, and many different choices of 
qt, as indicated. In each case, two straight lines are shown: the result of des Cloizeaux Eq. {TJ, qLS(q) = %+ (2/3)(q£ p )~ 1 , 
and empirical fits, qLS(q) — ty + 2.35(q£ p )~ 1 for d = 2 and qLS(q) = ty + 1.9(q£ p )~ 1 for d = 3, respectively. 



Another interesting behavior is the apparent power 
law, {R 2 g± ) oc (fe p /k B T)~ 2 - 5 , seen for f£ p /k B T > 1. 
However, we warn the reader to take this seriously: a 
closer look reveals a slight but systematic curvature, and 
a plot versus flp/ksT on a linear rather than a logarith- 
mic scale reveals that this apparent power law is noth- 
ing but the onset of an exponential decay (Figs. [SJ:,d): 
indeed, already from the partition function, Eq. (|62[) 
we recognize that for large forces the chains will be 
stretched out almost completely like rigid rods, and the 
few remaining kinks are suppressed exponentially when 
f£ p /k B T^l. 

For completeness, we show the corresponding simula- 
tion data for the relative extension (X}/L vs. ftp/ksT 
in Fig. [9] (related more extensive data for other values 
of L and qb can be found in Ref. and in Fig. \W\ 

we present the corresponding data for the longitudinal 
component (R 2 g 1 1 ) /L 2 of the gyration radius in the di- 



rection of the force. While for very small forces one ex- 
pects a nonzero plateau (unlike (X)/L, which vanishes 
as / — > 0), corresponding to the gyration radius square 
component of an unstretched chain, for large / another 
plateau means that the chain has been stretched out fully 
to a rod of length L. In between these two plateaus, the 
Pincus blob behavior is seen rather clearly, for the flexible 
chains. 

Fig. 11 now shows typical data for S±(q±) vs. q± and 
Fig. 12 the corresponding data for S^^h) vs. qii, focus- 
ing again on those selected values of qb that were used 
already in Figs. 8-10. As expected from Eq. (|4"T)) . the per- 
pendicular structure factor is similar to the case without 
stretching force; the plateau at small q, where S±(q±) 
deviates only very little from unity, gets more extended 
with increasing /, reflecting the decrease of (R 2 j_) with 
/ (Fig. 8). This decrease, of course, is more pronounced 
for stiff chains than for flexible chains at the same value 



13 




FIG. 8: Log-log plot of (R^±) / '(R%,±)o vs. f£ p /k B T, for d = 2 (a) and for d = 3 (b), using L = N£ b = 25600 in both cases, for 
two choices of qt- In parts (c,d) the same data as in (a,b) for large flp/ksT are replotted choosing a linear scale for flp/ksT, 
to show that the behavior is compatible with an exponential decay in d = 3. 



of /, since we have seen (see Fig. 8 and the discussions in 
Ref. [l3j]) that the proper control variable is not f/kgT 
but rather fip/ksT. A further remarkable feature is 
the fact, that the power law-like decay of S±(q±) with 
q±, that for flexible chains can be observed until S±(q±) 
has decayed up to S±(q±) ~ 1CP 4 as b < 1.5, for stiff 
chains extends only to S±(q±) ~ 10~ 3 in d = 3 and 
S±(q±) « 10~ 2 in d = 2 respectively. As expected, the 
theory of Benoit et al. [4(| which extended the description 
of scattering from Gaussian chains to elastic stretching 
deformations, can only be applied if q±£ p <IC 1, and when 
(Rg j_) exceeds £p only by few orders of magnitudes, the 
applicability of Eq. ([4"T| is correspondingly restricted. In 
fact, noting that, for = 0.05, l v jlb ~ 5.9, we conclude 
that b = 1.5 means f£ p /ksT w 2.4, and Fig.[8]shows that 
in this case indeed (i? 2 ± ) is about an order of magnitude 
smaller than for / = 0. Assuming that one can represent 
a stiff chain as a sequence of rods of length l v such that 
n p £p = Nib, and stating in the spirit of Eq. ([To]) that at 
large q± interference effects of different rods can be ne- 
glected, one would expect that for q±£ p ~ 1 one obtains 
a scattering of the order of S±(q±) ~ n~ 1 , independent 
of q±. This (admittedly rough) argument would qualita- 
tively explain the systematic increase of the plateau S(q) 
in Fig. 5 and S±(q±) in Fig. 11 with increasing chain 



stiffness. 

Even more interesting is the behavior of Si i (q> i ), 
Fig. 12. The rapid increase of (i? 2 ||) with increasing 
stretching force has the consequence that S , ||(g[|) devi- 
ates from unity for smaller and smaller qi | . While for 
small / just a shoulder develops, before (at large qn) the 
behavior is similar to that of S±(q±), for large / pro- 
nounced oscillations develop. As pointed out already by 
Pierleoni et al. [44[ for the case of fully flexible chains 
under stretch, this behavior can be attributed to the 
fact that the chain behaves like an elastically stretched 
string. We shall discuss this behavior in more detail be- 
low. Here we only note that the maxima of these oscil- 
lations decay according to a power law, which is similar 
to the power law of S±(q±) in the intermediate range 
of q±. The minima of Si^gii), as well as SW^ii) itself 
at larger values of q\\ where the oscillations of q\\ have 
decayed, show a slow further decrease with qn. While 
for very large qu = q± but not very strong stretching 
(b < 1.5) we have S\\(q\\) = S±(q±), if the chains are 
flexible (qt, = 1, qj, = 0.4), this is not the case for stiff 
chains: Si^gii) <C S±(q±) for qu = q± then. 

In order to understand these results more quantita- 
tively, we first report S±(q±) in form q±S±(q±) versus 
q± and compare to the Debye function, Eq. (|47[) , but us- 



14 




FIG. 9: Log-log plot of the relative extension {X)/L versus the scaled force flp/ksT for L = 25600 in d = 2 (a) and d = 3, for 
the same choices of qb as shown in Fig. [8j The straight lines indicate the linear response regime and the Pincus-blob regime, 
respectively {Eq. (|59[) }. In neither case has the Kratky-Porod model (curves labeled as K-P-model, cf. Eqs. (|54|l . (|55[) ) a well- 
defined regime of validity (this is expected since for large flp/ksT the approach to saturation, (X)/L — > 1, is exponentially 
fast rather than proportional to (fcsT/ fip) 1 ^ 2 , cf. Fig. [8}. 
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ing the value (R 2 g j_) as observed in the simulation (rather 
than any theoretical prediction for it). Fig.fTBlshows that 
the Debye function works surprisingly well: for d = 3 the 
slope 2 — 1/v indicating non-Gaussian behavior is seen 
for qh = 0.4 only for weak stretching (6 = 1.001 and 
1.003), while for larger stretching forces a horizontal part 
in the plot q 2 j_S±(q) has developed. Also in d = 2 the ex- 
cluded volume regime, where the slope 2 — 1/v is compat- 
ible with the data, is pronounced only for rather flexible 
chains (such as qt, = 0.4, Fig. [TBI:) while for stiff chains in 
d = 2 (such as qb = 0.01) excluded volume effects show 
up in S±(q±) only for extremely weak stretching (such 
as b = 1.0001, i.e. //fcsT = 10~ 4 ). for stronger stretch- 
ing of semiflexible chains in d = 2 the Debye function 
seems to describe the data for small q± (q± < 10~ 2 ), but 
then a crossover to a behavior Sj_(qj_) ~ const and hence 



q±S±(q±_) oc q\ sets in. 

The behavior of Si \{q\\) when plotted in the form 
q^S\\(q) vs. <7|| is particularly striking (Fig. 14). Again 
excluded volume effects (described in this representation 
by a slope 2 — 1/v again) are pronounced only for very 
small forces, while for somewhat larger forces (such as in 
d = 3 for b > 1.070 for q b = 0.4 and for b > 1.015 for 
qb = 0.05, and in d = 2 for b > 1.015 for q b = 0.4 and 
for b > 1.005 for qb = 0.01), the oscillatory behavior of 
the structure factor, as described by the Debye function 
with complex X\\ {Eq. (|46p }. sets in. In order to inter- 
pret this behavior in more detail, we write X\ i = a + ic, 
where a = q^((X 2 } — (X} 2 )/2 and c = q\\(X), to rewrite 
Eq. (g7|) as follows 
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FIG. 11: Log-log plot of the structure factor S_i_(q±_) vs. gx for d — 3 and two choices of the stiffness parameter qb, qt = 0.4 (a) 
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b = exp(//fc_gT) are included, as indicated. Chain length is L — 25600 throughout. 



S||(<7||)=2 



exp(— a) [(a 2 — c 2 )cosc— 2acsinc] 



(a 2 +c 2 ) 2 
I 



(63) 



We recognize that there are two rather distinct parts, 
an exponentially damped oscillatory part and a "back- 
ground part" which survives when the oscillatory part 
has died out. For a ^> 1 this background part can be 
written as (note that for large stretching there is a regime 
where a 2 < c 2 ) 



- 1 



(64) 



In the regime where the oscillations have died out again 
and there is hence a fiat part of Sm^m), independent of 
<7|| again, the structure factor hence measures the rela- 
tive fluctuation in the length of the strongly stretched 
polymer. 

Let us now consider the oscillatory part of Eq. (1631 . 
Since we are in a regime where a 2 <C c 2 , the maxima are 
reached when cos c = — 1 and in the regime where a«l 



1 we hence have 

4 _ 4 
~ q^X) 2 ' 

(2m + 1)tt, m = 0, 1, . . 



(65) 



and hence cxp(— a) 

q\\(x) 

and hence q 2 S™ aK ps const, as observed from the full 
calculation of Eq. (|63|) . and the simulation. 

When we compare these results to the scattering from 
the rigid rods, however, Eq. (|30p predicts maxima that 
are undamped and minima that are strictly zero, so for 
increasing qu the oscillations continue forever. However, 
this is a result for a rod that has a strictly fixed length 
L m d, while the polymer under strong stretch (with ex- 
tension (X) such that 1 — (X) / L <C 1) still is only similar 
to a rod of fluctuating length, and this in fact is borne 
out by the structure factor at large q {Eq. 

In any case the success of the Debye function, 
Eqs. (Unj), (jSSJ) for the description of the scattering from 



16 



(a) 



(b) 




10" 



S 10" 2 



10"' 



10- 



10 




L = 25600 
b 

1 .000 
1.005 
1.025 

1.070 o\*SM 
1.150 *i i 
1.250 , □ 



(c) 




1 71 



111 



(d) 



CO 




111 



FIG. 12: Log-log plot of the structure factor Sii^ii) vs. qu for d = 3 and two choices of the stiffness parameter qb, qb ~ 0.4 (a) 
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strongly stretched chains in both d = 2 and d = 3 dimen- 
sions is very remarkable, since it is derived from Gaussian 
chain statistics |40j . and we have seen that in d = 2 in 
the absence of stretching forces Gaussian chain statistics 
does not work in d = 2, irrespective of chain stiffness. 



At the end of this section, we emphasize that the ex- 
amples given for the success of the Debye function for 
stretched chains, as derived by Benoit et al. (40[, are not 
accidental, but this behavior is typical for a wide range 
of chain stiffnesses. As an example, we show further 
data for S±(q±) in both d = 2 and d = 3 and various 
other choices of the stiffness parameter qi, in Figs. 15, 
16. Whenever the plots indicate a well-defined plateau, 
one can extract an estimate of (R 2 j_) from it (note that 
the actual values of (Rg ±) that were independently esti- 
mated were used to predict the Debye functions as shown 
in Figs. 15, 16.). As has been shown already in Fig. 8, 
the radii (R 2 ± ) do have a broad regime of forces f/ksT 
where excluded volume effects ("Pincus blob'" - behav- 
ior) prevail, so the success of the Debye function must 
not be over-emphasized, it does not mean that the chain 
conformation follow Gaussian statistics. 



VI. CONCLUSIONS 

In this paper we have presented a comparative simu- 
lation study of the single-chain structure factor S(q) for 
variable stiffness of the macromolecules in both d = 2 
and d = 3 dimensions, both for coils in equilibrium in 
dilute solution under good solvent conditions, and for 
polymers under the influence of a stretching force. Char- 
acteristic linear dimensions of the macromolecules that 
are needed in the theoretical interpretation of S(q), have 
in our Monte Carlo simulation always been estimated di- 
rectly and hence independently, such as the mean square 
gyration radius (Rg) and the persistence length £ p . In the 
presence of stretching forces, the extension (X) in the di- 
rection of the force (as well as fluctuations (X 2 ) — (X) 2 , 
and components of the gyration radius (Rg±) and (R^) 
have been obtained as well. Of course, as usual the simu- 
lations are performed for the strictly monodisperse case, 
the number of bonds N and hence also the contour length 
L = N£{, of the chain molecules are known input param- 
eters of the simulation. In this respect, a more definite 
interpretation of the outcome of the simulations can be 
expected, than would be expected for corresponding ex- 
periments (where polydispersity is a problem, and of- 
ten the average contour length is not a priori known but 
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FIG. 13: Log- log plot of q±S±(q±) vs. q± for d = 3 and two choices of the stiffness parameter (ft, (ft = 0.4 (a) and qt = 0.05 (b), 
and also for d — 2 and two choices of Qb, namely qt = 0.4 (c) and qt = 0.01 (d). The curves are the Debye function, Eq. (|47[) 
with X± — |g5_(7?g ± ) in d = 3, and X± = 3q±(R^ x) in d = 2. In each case several choices of 6 = exp(//fc_eT) are included. 
Chain length is L = 25600 throughout. 



must also be extracted from fitting suitable experimen- 
tal data). Of course, the drawback of our Monte Carlo 
simulations on a lattice is the highly idealized charac- 
ter of such coarse-grained model as used here, the self- 
avoiding walk with additional energy penalty for making 
kinks. Nevertheless, the comparison between our simu- 
lation results for the mean square gyration radius versus 
the number of "Kuhn segments" uk with corresponding 
experimental data (Fig. 1) is very encouraging: one not 
only notes a striking similarity between simulation and 
experiment, but we stress that also the same range of di- 
mensionless variables (uk and (R 2 g ) / (2£ p ) 2 ) is accessible. 
The simulation has the bonus that directly single-chain 
properties are accessible (no extrapolation as a function 
of the concentration c of the solution towards c — > is 
required and by changing the energy parameter et/ksT, 
that describes the cost of making a kink along the walk, 
in units of the thermal energy, the stiffness is easily con- 
trolled. In experiment, stiffness can only be widely varied 
by combining data for polymers with different chemical 
structure. 

From our data we have confirmed a conclusion drawn 
already from our previous study of mean square end-to- 
end distances, namely that in d — 2 a direct crossover 



occurs from rod-like behavior to self-avoiding walks, with 

a scaling (Rg) oc £p^ 2 L 3 ^ 2 , without the existence of any 
intermediate regime with Gaussian behavior (Fig. 3). In 
d = 3, however, such an intermediate regime has been 
found, Fig. 2, for 1 < n K < n* K oc (£ P /D) C , where D 
is the local chain diameter and the exponent £ is in the 
range 1.5 < £ < 2. Thus, there is no universal value 
where excluded volume effects set in, but rather n* K — > oo 
for £ P /D — !• oo. 

In the equilibrium structure factor S(q), in the absence 
of stretching forces, correspondingly several regimes can 
be distinguished. For small enough q, the standard 
Guinier behavior always occurs, which contains the in- 
formation on (Rg), of course. For d = 2, one then always 
has the excluded volume regime (for long enough chains), 
S(q) oc q -4 / 3 , and possibly (for rather stiff chains) a 
crossover to rod-like behavior (S(q) oc g _1 ) sets in grad- 
ually. A Gaussian behavior S(q) oc q~ 2 is never seen, 
unlike the case d — 3, where this behavior does become 
visible for very stiff chains (before for still larger q the 
rod- like behavior starts). The excluded volume power 
law, S(q) oc q" 1 ^ with v sa 0.588, is only visible for not 
very stiff chains (if chain lengths L < 25600 are analyzed, 
as done here: if L — > oo, this power law would emerge for 
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FIG. 14: Log-log plot of q?\S\\ (q\\) vs. for d — 3 and two choices of the stiffness parameter qb, namely qt = 0.4 (a) and 
qs = 0.05 (b), and also for d = 2 and two choices of qb = 0.4 (c) and qb = 0.01 (d). The curves are the Debye function, 
Eq. (|63p with complex X\\ = q^((X 2 ) — (X} 2 )/2 + iq\\ < X >, and the straight line shows the excluded volume power law 
(slope — 2 — In each case several choices of b = exp(//fesT) are included. Chain length is L = 25600 throughout. 



any finite value of the persistence length, of course). This 
pattern of behavior (Fig. 4), of course, could have been a 
priori expected, but we are also able to show via Kratky 
plots (qLS(q) vs. qL) that for semiflexible chains in d = 3 
the expressions derived by Kholodcnko and by Stepanow 
provide a quantitatively accurate description. It is found 
that for large q this quantity qLS(q) settles down at n, 
unlike the behavior predicted for flexible chains (the De- 
bye function predicts qLS(q) oc q^ 1 for large q). How- 
ever, the onset of the plateau occurs gradually in the 
decade 1 < q£ p < 10; thus the onset of the plateau allows 
an estimation of £ p only somewhat roughly. The peak 
position of the Kratky plot (Figs. 5,6) reflects the theo- 
retically expected scaling of the gyration radius with L 
and l p , even though in the Kratky plot (Fig. 5) direct 
evidence for excluded volume effects seem to be minor. 

Des Cloizeaux [23| derived from the Kratky-Porod 
model that for L —> oo one should have LqS(q) = 
7r + const(g€ p ) _1 , with the constant being predicted to 
be 2/3 {Eq. ([TJ}. Unfortunately, this result is at vari- 
ance with our numerical results (Fig. 7). The reason for 
this problem is still not clear. 

Turning to the behavior of chains under the influence 
of stretching forces, we have shown that for weak forces, 



where linear response holds, excluded volume effects in- 
validate the Kratky-Porod model completely in d = 2 
dimensions, and one typically observes a broad range 
of forces where the extension versus force relation is a 
power law, and also (Rg±), (X 2 ) — (X) 2 are found to 
scale like (/ip/hsT) 1 ^- 2 in this "Pincus blob" regime. 
In d = 3, dimensions, however, a Pincus blob regime also 
exists, but its observability also is restricted. For large 
ksT ' I f ', however, the (continuum) Kratky-Porod descrip- 
tions is not valid for our discrete lattice model either: it 
is found that then (R^±) and (X 2 ) — (X) 2 decrease like 
exp(-//fc B T) for f/k B T » 1. 

Although the excluded volume effects show up rather 
clearly in the chain extensions and gyration radii compo- 
nents, it turns out that Benoit's extension of the Debye 
formula to stretched chains [4(J is a surprisingly accurate 
description of both the transverse (S±(q±)) and parallel 
parts (5||(g||)) of the structure factor. The oscillatory 
behavior of <S'||(q , ||) for strongly stretched chains shows 
that their conformations resemble a string of elastically 
coupled particles. Thus, if measurable, the structure fac- 
tor of stretched chains would add valuable details to the 
picture of their conformations. 

As we have emphasized in our paper, the statistical 
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FIG. 15: Log-log plot of q^S±(q±) vs. qx in d = 3 dimensions for qt — 0.2 (a), 0.1 (b), 0.03 (c) and 0.01 (d), for several choices 
of the force parameter b — exp(/ /ksT), as indicated. The curves are the Debye function, Eq. (|47|l with X± = ^q±{Rg ±). 
When the data settle down at a horizontal plateau, it yields an estimate of 4/(3(JJg 



mechanics of semiflexible polymers has been a longstand- 
ing and controversial problem of polymer science. The 
subject is of great relevance for biopolymers, but also of 
broad interest in material science. We expect that the 
present study will be useful both for the interpretation 
of experiments and stimulate further theoretical studies, 
such as of the interplay between solvent quality and chain 
stiffness. 
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Appendix A: Scattering function of random walk 
chains under constant pulling force 

The paper by Benoit et al. [io| uses the distribution 
function of the end-to-end vector of a Gaussian chain of 
| i — j | repeat units for the calculation of the single chain 
scattering function. The result is the well-known Debye 
function. 

Following an idea from Doi's book [74| one can de- 
rive the diffusion equation yielding the end-to-end vector 
distribution in the following way. Assume a step-wise 
Markov growth of the chain 

P(R,N) = Y,P(R-h,N -l)p(h). (Al) 

i 

The sum over "i" goes over an isotropic bond vector set, 
i.e., both li and — Zj are members of the set. 

Expanding the right side to first order in N and to 
second order in R yields 
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Averaging the derivatives with respect to R a with a sym- 
metric bond probability p(U) (in the simplest case this is 
just one over the number of bonds) gives zero for the first 
derivative and 8 a pl 2 /3 for the second term, resulting in 
the diffusion equation 



dP 



G 



-V 2 P 



(A3) 



Solving this for a bulk chain gives the well known result 
for the end-to-end vector distribution 



P(R,N) = (2irNl 2 /3y 3/2 



exp 



3R 2 
~2Nf 



(A4) 



The Debye function can than be derived by averaging 
with this probability, as is done in the Benoit et al. paper. 

The derivation above is useful as a starting point for 
a calculation of the scattering function for a chain that 



is pulled. In this case the bond probabilities are not 
symmetric. For our model we have 
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(A5) 



(A6) 



for moves perpendicular to the pulling direction, in +X 
direction and in —X direction, respectively, where b — 
exp(flfkBT) with I — 1 is used as in the main text. 
When we now perform the expansion of Eq. (|A2[) and 
perform the average over the bond probabilities we obtain 
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For p + = P- = Po this reduces to the normal diffusion 
equation. This equation has to be solved with the bound- 
ary conditions 



P{R,0) = 8(R) 
P(R,N) -> for 



(A8) 



R 



Let us define D± = 2p l 2 , D\\ = {p + + p^)l 2 , and v — 
(p + — p-)l, so we have 

dP _ BP 1 fd 2 P d 2 P\ 1 d 2 P 
M~~ V dX + 2 - 1 [dY 2 + dZ 2 J + 2 ^dX 2 ' 

(A9) 

These are three diffusion processes in the three Carte- 
sian directions, X is parallel to the force, Y and Z are 
perpendicular. The solutions for the perpendicular di- 
rections are the same as for the force-free case. For the 
parallel direction we have an altered diffusion coefficient 
and a drift part to the process, i.e., a Gaussian diffusion 
around a deterministic drift. The complete solution to 
Eq. (|A9|) is therefore given by 



P(X, Y, Z, N) 



1 



1 



2nND± v /27r7VZ>| i 



(A10) 

For po = p+ = p~ = 1/6 we obtain back the force free 
solution. To calculate the scattering function we follow 
the procedure employed in the calculation of the Debye 
function in the force free case. 



S(q) 



1 



N 2 ^ 



(All) 



can be calculated when we assume a continuous chain 
model (i.e. only look at distances much larger than the 
lattice constant) so that the distribution for the is 
given by the Gaussian distribution we just calculated. 

s ® = ^ E / <? F n p (?ii > i * - 3 \y^ 3 ( Ai2 ) 



where we have P(rij, | i — j |) = Py (Y, \ i — j \)Pz{Z, | 
i — j \)Px(X, | i — j |) and Py and Pz have the same 
functional form and all Pi are normalized to one individ- 
ually. 



where p^j = Yey + Z&z- So we have to evaluate 
S(q Y ,q z ) 

If 1 Y 2 



N 2 ^J ~ ^\i- 3 \D A 



dZ- 



V2tt 1 * — j 1 D 



_g 3|i-3|-D_L g'9 z 



resulting in 



1 ^ _ g j^ k iwi 



(A14) 



(A15) 



This is evaluated by a continuum approximation for the 
two sums, — ¥ du and y\ — > dt> which finally 
yields: 



+ 



8 



Nq 2 ,D ± N 2 q 4 ,D 



e~^>— - 1 . (A16) 



Scattering in the parallel direction 



S (V\\) = I dX 



'■3 



which now results in 



2|i-j]-D M e n\\X 
(A17) 



e ^ cos(g||W | i - j |) . (A18) 



'■j 



Performing the final calculation again in the continuum 
approximation gives 



4 



N q*D 2 + 4v 2 q 2 
8 qf^-^q 2 



N 2 (qf,D 2 +4v 2 q 2 ) 2 



' e 2 cos(Nvq\\) - 1 



32 9m-^II w 

^ PF^IF 6 "^ sin(iYw9||) (A19) 



Scattering in the perpendicular direction 

s (9±) = ^2 T.J dY J dz ( A13 ) 



This result determines our scattering functions with pa- 
rameters depending on the applied force /, so these equa- 
tions contain no free parameters. Note that both func- 
tions reduce to the Debye function for the force free 
isotropic case (v = 0, D± = Du = l 2 /3) as it should 
be, because the scattering function does not depend on 
the direction of the scattering vector in this case. 



22 



For fitting purposes it might yield better results to re- 
place some of the quantities by average values determined 
in the simulation. For the Y and the Z components sim- 
ple Gaussian statistics holds 



on 

(Y 2 ) 



(R 2 



9,YI 



(Z) 

(z 2 )-- 







NDj_ 
_ ND ± 
~ 6 



(A20) 
(A21) 

(A22) 



but in the force direction we can calculate from Eq. (|AlOj) 
for the moments of the end-to-end distance 



(X) = Nv 
(X 2 ) = 
(AX 2 ) = 



N 2 v 2 



(R 2 



g,xi 



ND\\ 
ND\ 
6 



■NDu 



D\\N 2 v 2 
12 



(A23) 
(A24) 
(A25) 

(A26) 



When we now rewrite the scattering functions in terms 
of moments of the end-to-end vector or gyration tensor, 
we obtain 



This gives the same formula as Eq. (|63j) when Xu in 
Eq. (gBJ is written by X\\ = a + ic with a = q?, (AX 2 ) /2, 
and c = q\\(X). For (X) = and (AA 2 ) = (X 2 ) = 
6(R 2 x ) = 2(i2g) this again reduces to the Debye func- 



tion, Eq. (|22j>. 



Scattering in d = 2 



Both scattering functions as calculated in Eqs. (|Al6l) 
and (|A19[) remain formally unchanged. However, the 
probabilities for the single steps change to 



Po 

P+ 
P- 



b 2 - 


f 26 H 


- 1 




b 2 




b 2 - 


f- 26H 


- 1 




1 




b 2 - 


F 26 H 


- 1 



(A29) 



S(qj 



(A27) 



e - 1 



This is the Debye function with the appropriate pref- 
actors, because for (R 2 ± ) = 2/3(R 2 ) it reduces to the 
standard function, Eq. (|2"2"j) . For the scattering parallel 
to the pulling direction we can write 



The parallel and perpendicular diffusion coefficients as 
well as the drift velocity still have the same functional 
dependence on these probabilities. However, introducing 
the chain extensions into Eqs. (|A16[) and (| A19I) for d = 2 
changes the prediction Eq. (|A28[) for the perpendicular 
scattering to 



5(5, | ) (A28) S{q±) 2 (A30) 



, g H (AX } 3ql (R 2 x ) + (3ql(R 2 ± )) 2 V 6 V 

qfAAX 2 ) 2 + Aq 2 (X) 2 ^9,^ \1X\g,±JJ 



qUAX 2 ) 2 -4q 2 (X)' / ^ 



q?,(AX 2 )(X) g f|(AX^> *<l±\X-g 

- 3 \qUAX 2 ) 2 +4q 2 (X) 2 ) 2 e ~ 2 ^iW)- mulaasEq.gTl. 



whereas it leaves Eq. (IA28[) unchanged. Taking X± = 
3q 2 [(R 2 ± ), the expression of Eq. (|A30[) has the same for- 
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